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ABSTRACT 

The distribution of galaxy properties in groups and clusters holds important in- 
formation on galaxy evolution and growth of structure in the Universe. While clusters 
have received appreciable attention in this regard, the role of groups as fundamental to 
formation of the present day galaxy population has remained relatively unaddressed. 
Here we present stellar ages, metallicities and a-element abundances derived using Lick 
indices for 67 spectroscopically confirmed members of the NGC 5044 galaxy group with 
the aim of shedding light on galaxy evolution in the context of the group environment. 

We find that galaxies in the NGC 5044 group show evidence for a strong rela- 
tionship between stellar mass and metallicity, consistent with their counterparts in 
both higher and lower mass groups and clusters. Galaxies show no clear trend of age 
or a-element abundance with mass, but these data form a tight sequence when fit 
simultaneously in age, metallicity and stellar mass. In the context of the group envi- 
ronment, our data support the tidal disruption of low-mass galaxies at small group- 
centric radii, as evident from an apparent lack of galaxies below ^ 10^ Mq within 
~ 100 kpc of the brightest group galaxy. Using a joint analysis of absorption- and 
emission-line metallicities, we are able to show that the star-forming galaxy popula- 
tion in the NGC 5044 group appears to require gas removal to explain the '^1.5dex 
offset between absorption- and emission-line metallicities observed in some cases. A 
comparison with other stellar population properties suggests that this gas removal is 
dominated by galaxy interactions with the hot intragroup medium. 

Key words: galaxies: clusters: NGC 5044 group - galaxies: distances and redshift - 
galaxies: fundamental parameters - galaxies: formation - galaxies: evolution 



1 INTRODUCTION 

Large photometric surveys have established that the ob- 
servable galaxy population forms two distinct sequences in 
colour-magnitude space: a red sequence of predominantly 
early-type galaxies and a star-forming, late-type dominated 
blue sequence (e.g. Hogg et al. 2002; Blanton et al. 2003; 
Baldry et al. 2006). Recent observational evidence suggests 
that growth of the red sequence since z ~ 1 is dominated 
by increasing numbers of low- luminosity (<L*) red galax- 
ies with cosmic time at the expense of the low-luminosity 
star-forming galaxy population (e.g. Bell et al. 2004; Bundy 
et al. 2006; Brown et al. 2007; Faber et al. 2007). Hence, our 
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understanding of growth and evolution in the red sequence 
is closely related to our knowledge of processes responsible 
for quenching star formation in blue-sequence galaxies, en- 
abling their subsequent evolution onto the red sequence. 

The dichotomy of observed galaxy colours carries with 
it an increasingly well-described environmental dependence 
whereby galaxies inhabiting high-density regions are on av- 
erage redder, older, more concentrated and more massive 
than galaxies found in the field (e.g. Oemler 1974; Dressier 
1980; Goto et al. 2003; Blanton 2005; Blanton & Berlind 
2007; O'Mill, Padilla & Lambas 2008; Cooper et al. 2008). 
This strong correlation of red-sequence galaxies with local 
galaxy density has led to a significant body of work devoted 
to exploring the environmental factors responsible for creat- 
ing these predominantly red populations from the blue, star- 
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forming population of field galaxies. Of the possible mecha- 
nisms for suppressing star formation, ram-pressure stripping 
(e.g. Gunn & Gott 1972; Hester 2006), mergers (e.g. Toomre 
& Toomre 1972; Negroponte & White 1983; Somerville & 
Primack 1999; Mcintosh et al. 2008), strangulation (e.g. 
Larson, Tinsley & Caldwell 1980; Balogh & Morris 2000), 
tidal interactions and harassment (e.g. Farouki & Shapiro 
1981; Moore et al. 1996) are the most frequently invoked 
to account for the increasing number density of red galax- 
ies. Identifying which of these possible processes represents 
the dominant evolutionary path, however, remains an active 
field of research. 

The Sloan Digital Sky Survey (SDSS) has proved funda- 
mental in exploring the environmental dependence of galaxy 
properties, and has been used to show that the transfor- 
mation of galaxies from blue to red takes place over rela- 
tively long timescales, > 1 Gyr, through the slow exhaustion 
of gas available for star formation, rather than rapid, vio- 
lent removal (e.g. Kauffmann et al. 2004; van den Bosch 
et al. 2008). The dominance of so called "strangulation" 
(Balogh & Morris 2000) precludes cluster-cores as domi- 
nant sites for blue galaxy transformation as disruption pro- 
cesses in high-density regions, e.g. ram-pressure stripping 
and harassment, act over relatively short timescales (e.g. 
Abadi, Moore & Bower 1999; Quilis, Moore & Bower 2000). 
This supports earlier work from the 2dF Galaxy Redshift 
Survey (2dFGRS) which found that the fraction of star 
forming galaxies is lower with respect to the field even at 
the relatively low surface densities of loose groups, ~1- 
2 galaxies/Mpc^ , i.e. well outside dense cluster cores (e.g. 
Lewis et al. 2002; Poggianti et al. 2006). 

While the majority of our understanding of the varia- 
tion of galaxy properties with environment hinges on large 
photometric surveys, an alternative approach is the use of 
spectral absorption features as detailed tracers of galaxy 
stellar populations. In many ways spectral analyses provide 
a vastly increased level of detail for determining paths of 
galaxy formation and evolution relative to photometry, but 
they are often hindered by smaller sample sizes and rela- 
tively high observational expense. Nevertheless, absorption- 
line studies have met with great success in tracing the stellar 
populations of individual galaxies (e.g. Proctor et al. 2004; 
Thomas et al. 2005; Sanchez-Blazquez et al. 2006; Peletier 
et al. 2007; Spolaor et al. 2008), clusters (e.g. Smith, Lucey 
& Hudson 2007; Smith et al. 2008a,b; Trager et al. 2008) 
and large survey samples (e.g. Gallazzi et al. 2006; Jimenez 
et al. 2007; Proctor et al. 2008). 

To date, spectroscopic studies with a particular focus 
on environment have primarily made use of cluster galaxy 
populations. Caldwell, Rose & Concannon (2003) used a 
sample of nearby early-type galaxies in the Virgo Clus- 
ter, as well as a field-galaxy sample, to probe galaxy stel- 
lar populations to low masses (stellar velocity dispersions, 
0", < 100 kms~^). Caldwell et al. (2003) found a greater in- 
trinsic scatter in the properties of low-mass galaxies, but 
made no particular effort to examine the spatial distribution 
of galaxy populations within the Virgo Cluster. In similar 
work. Smith et al. (2007) and Smith et al. (2008a,b) have 
examined the stellar populations of dwarf galaxies in the 
Coma Cluster and Shapley supercluster, in both instances 
finding that the scaling relations of low-mass populations 
are generally consistent with the relations of higher-mass 



E and SO galaxies, although with the same increased scat- 
ter noted by Caldwell et al. (2003). In the case of Coma, 
Smith et al. (2008a) find a radial dependence of both age 
and metallicity, where galaxies at larger projected radii are 
both younger and higher metallicity, suggestive of recently 
quenched star formation among galaxies entering the Coma 
cluster outskirts and consistent with earlier spectral studies 
of Coma cluster galaxies (Caldwell et al. 1993; Mobasher 
et al. 2001; Carter et al. 2002). 

These studies, however, only serve to probe rare high- 
density cluster environments, neglecting the much more 
common group environment. In an attempt to fill this gap, 
here we undertake a detailed spectroscopic study of the 
galaxy population in the NGC 5044 group with the aim of 
describing galaxy evolutionary histories and reconstructing 
the formation history of the group as a whole. In a previ- 
ous paper (Mendel et al. 2008; hereafter Paper I), we used 
new, deep spectroscopic observations in conjunction with 
archival velocity information to more than triple the num- 
ber of velocity-confirmed group members. This greatly im- 
proved sample has allowed us to characterise the dynamical 
state of the group, which we find to be relaxed from a com- 
bination of X-ray and dynamical indicators. Here, we turn 
our focus to analysing the stellar populations of galaxies in 
this redefined group with a focus on establishing the details 
of galaxy evolution in the context of the group environment. 

The structure of this paper is as follows: Section 
[2] describes the spectroscopic data used from both new 
AAOmega observations and the 6dF Galaxy Survey. Here 
we also describe the method used to derive kinematic prop- 
erties. In Section[3]we discuss the measurement of Lick index 
absorption features and their interpretation in terms of age, 
metallicity and Q-element abundance (Section |4|. In subse- 
quent sections we discuss the stellar population properties 
of NGC 5044 group galaxies both as an independent galaxy 
sample (Section [5}, and in the context of their place in the 
group as a whole (Section [6)l . In Section [7] we describe the 
distribution of emission-line properties for NGC 5044 group 
galaxies. In Section [8] we discuss our stellar population find- 
ings in relation to the dynamical description of the group 
presented in Paper I, as well as their relevance to the for- 
mation and evolution of the NGC 5044 group and galaxies 
therein. 

Throughout this paper we assume Ho = 
70kms~^Mpc~^ where relevant. We adopt a distance 
modulus for the NGC 5044 group of (m - Af)o = 32.31 
(28.99 Mpc) measured using surface brightness fluctuations 
(Tonry et al. 2001) with corrections applied to adjust for 
the improved Cepheid distance measurements of Jensen 
et al. (2003). Where specified, correlations have been 
measured using the Kendall rank-order correlation test 
with significance quoted from the two-tailed p-value. 



2 DATA 

Here we use the NGC 5044 group sample defined in Paper 
I, which consists of 111 spectroscopically-confirmed group 
members compiled using newly obtained redshifts from 
AAOmega supplemented using the 6dF Galaxy Survey Data 
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Release 2 (6dFGS DrS Jones et al. 2005), HI observations 
(Kilborn et al. in preparation) and the Nasa/IPAC Extra- 
galactic Database (NED) sources. Of these 111 galaxies, 95 
have spectra available from AAOmega and 6dFGS observa- 
tions which can be used to measure absorption indices (71 
from our new AAOmega observations and 24 from 6dFGS 
DR2). These spectral data are summarised below. 



2.1 AAOmega spectroscopy 

New observations of NGC 5044 group galaxies were carried 
out using the AAOmega multi-object spectrograph on the 
3.9m Anglo- Australian Telescope (AAT) in Siding Spring, 
Australia. Target galaxies were taken from the list of po- 
tential NGC 5044 group members of Ferguson & Sandage 
(1990). Medium resolution 580V and 385R gratings were 
used, yielding FWHM dispersions of 3.5 A and 5.3 A for 
blue (AA3700-5700A) and red (AA5700-8800A) spectra re- 
spectively. 

Galaxies were separated into high- and low-luminosity 
subsamples using an apparent magnitude cut of B = 15 
and observed separately in order to limit the effects of scat- 
tered light from the brightest sources. Observations were 
carried out 30 minutes at a time with a total of 3 to 4 ob- 
servations per plate configuration. In addition to fibres as- 
signed to galaxies, in each configuration an additional 30 to 
35 "empty" fibres were assigned to be used for sky subtrac- 
tion. There are several relatively wide defects in the blue 
arm CCD of AAOmega, and to limit the effects of these 
on our final coadded science spectra we offset the central 
wavelength on each night of observations by 100 A, ensur- 
ing that usable data were obtained from at least two out of 
three nights at every wavelength. Our total integration time 
for high- and low- luminosity galaxies was 4 and 17 hours 
respectively. 

In addition to our primary (galaxy) targets we have 
observed a set of calibration stars. These include stars com- 
mon to the Lick stellar library (see i]3.2p and a set of spec- 
trophotometric calibrators which we can use to roughly flux 
calibrate the AAOmega data. 

Data were reduced using a combination of iraf rou- 
tines and the 2dfdr reduction pipeline supplied and main- 
tained by the Anglo- Australian Observatory. Prior to re- 
duction with 2DFDR an average bias was subtracted from all 
frames and bad pixel columns were identified and repaired. 
We then used 2dfdr to reject cosmic rays, identify fibre 
apertures, flat-field, throughput calibrate and sky subtract 
the data. The final coaddition of science frames was carried 
out after excluding any single exposures with very low flux 
relative to the other frames. 

In this work, we consider only those galaxies identifled 
as group members in Paper I, and refer the reader to that 
work for the details of how group membership was deflned. 
We exclude from subsequent analyses any galaxies with a 
signal-to-noise (S/N) of less than 12, as below this we con- 
sider stellar population measurements to be highly unreli- 
able. This cut excludes 19 galaxies, leaving a total of 52 from 
the original AAOmega sample described in Paper I. 
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2.2 6dFGS Spectral Data 

The 6dFGS covers 17046 deg^ of the Southern sky using 
the Six-Degree Field (6dF) multi-flbre spectrograph on the 
UK Schmidt Telescope. The survey targets are selected from 
the 2MASS Extended Source Catalogue (2MASS XSC; Jar- 
rett et al. 2000), 2MASS and SuperCOSMOS to include 
all galaxies brighter than Ktot ~ 12.75 mag with wave- 
length coverage from AA4000-8400A. In Paper I we include 
16 unique 6dFGS galaxies and 8 galaxies that overlap with 
our AAOmega sample, all of which have spectra available. 
However, applying a similar cut in S/N as the AAOmega 
data above excludes an additional 9 of the 6dFGS galaxies, 
leaving a total of 15 for stellar populations analysis. 

The combined AAOmega and 6dFGS sample then con- 
sists of 67 confirmed group-member galaxies with S/N > 12, 
which we use for the remainder of this analysis. 



2.3 Recession velocity and velocity dispersion 
measurements 

Recession velocities and velocity dispersions were fit simul- 
taneously using the penalised pixel-fitting (pPXF) software 
of Cappellari & Emsellem (2004). Pixel fitting methods are 
notoriously sensitive to the templates used, so it is impor- 
tant to minimise the effects of so-called template mismatch 
on kinematic measurements. The pPXF code is an exten- 
sion of the Gauss-Hermite fitting method of van der Marel 
& Franx (1993) which allows for fitting to be carried out 
using a linear combination of template spectra, all but elim- 
inating template-mismatch-induced errors. 

As our input templates we use 50 stellar spectra 
from the MILES standard star library (Sanchez-Blazquez 
et al. 2006) spanning a broad range of spectral types. The 
particular choice of stars used to construct templates has 
only a minimal effect on the results derived using pPXF as 
long as they span a similar range of spectral types. We have 
investigated this briefly by fitting a sub-sample of our galaxy 
data with sets of 50 randomly chosen stars from the MILES 
library and find typical deviations of ~4-5 kms~^in velocity 
dispersion. 

In order for pPXF to provide meaningful velocity dis- 
persion outputs we need to account for the difference in in- 
strumental dispersion between our AAOmega spectra and 
the MILES stellar library. Overlap between the Lick and 
MILES stellar libraries is significant, so we have used the 
Lick standard stars observed throughout our AAOmega ob- 
serving run to effectively measure the wavelength depen- 
dent broadening offset between AAOmega and the MILES 
library. In total we use 5 Lick/MILES stars to measure the 
median broadening difference. We find that scatter in our 
broadening measurements is significant for wavelengths be- 
low 4700 A, but above this wavelength the dispersion offset 
is relatively stable at ~57.7kms"^(2.3 A FWHM). Adding 
in quadrature the spectral resolution of the MILES library, 
2.3 A, gives an average spectral resolution for AAOmega be- 
tween 4700 A and 5400 A of 3.2 A FWHM. This derived res- 
olution is consistent with that found by Smith et al. (2007) 
for spectra taken using the 580V grating on AAOmega. 

Standard star observations for the 16 6dFGS galax- 
ies that we define as group members are unavailable, but 
we can measure the broadening offset between AAOmega 



4 Mendel et al. 



and 6dF using the 8 galaxies present in both sam- 
ples. In the AA4700-5400 A wavelength range we find a 
mean offset between AAOmega and 6dFGS spectra of 
103.6kms"^(4.1 A FWHM). Adding in quadrature the pre- 
viously measured offset between AAOmega and the MILES 
library in this spectral region we adopt a total broadening of 
118.6kms"^(4.7 A FWHM). The difference in aperture size 
between the 6dF and AAOmega fibres (6.7" vs. 2.1") will re- 
sult in an additional broadening difference between the two 
galaxy spectra, independent of instrument resolution, of or- 
der 5 percent (j0rgensen, Franx & Kjargaard 1995). We do 
not take this aperture effect into account when calculating 
our broadening as the statistical nature of the j0rgensen 
et al. (1995) correction will likely only serve to introduce 
scatter in our relatively small number of galaxies, particu- 
larly given our broad range of galaxy types. 

Prior to measuring kinematics using pPXF, MILES 
template spectra are smoothed to the AAOmega and 6dFGS 
resolution using the measured broadening offsets. The pPXF 
code allows for a penalty to be applied to the higher order 
terms of the Gauss-Hermite polynomials and here we adopt 
a moderately high bias of 0.8 throughout fitting in order to 
limit spurious results from undersampled or low S/N data. 
Random errors on our kinematic measurements are then cal- 
culated using a series of Monte Carlo simulations spanning 
a range of velocity dispersions and S/N. The scatter in our 
determination of instrumental broadening discussed above 
adds an additional, systematic uncertainty to our veloc- 
ity dispersion measurements of 22kms~^and 41kms~^for 
AAOmega and 6dFGS galaxies respectively, which we in- 
clude in quadrature. 



2.4 Emission-line measurements 

We identify emission-line galaxies through examination of 
residuals to the best fit pPXF templates, and for 24 galaxies 
in our sample we find evidence for emission in some combina- 
tion of H^, [Om]A4959, [Om]A5007, [NI]a5198 and [NI]a5200. 

While we can estimate emission-line fluxes from the residu- 
als of our template fits, the underlying absorption in these 
regions is generally poorly constrained as emission contam- 
inated lines are masked in the template fitting process. 

A far more robust method is to fit for absorption- and 
emission-line kinematics simultaneously, which we carry out 
using the GANDALF software of Sarzi et al. (2006). This is an 
extension of the pPXF pixel fitting routine described in i]2.3l 
which fits simultaneously a set of absorption and emission 
templates, eliminating errors introduced by fitting emission 
line residuals separately. Fits to emission lines are conducted 
using a set of independent Gaussian profiles whose kinemat- 
ics can be varied as required. Here we treat the kinematics of 
nebular and Balmer-line emission separately, tying each to 
the kinematics of the dominant features in appropriate re- 
gions of our spectra (H/3 or Ha and [OIII]a5007 or [NII]a6583 
for Balmer and nebular lines respectively). Our choice to 
fit separately the kinematics of Balmer and nebular emis- 
sion lines has a negligible effect on our results as we find 
no galaxies with a significant kinematic difference between 
the two (where both are detected). There is generally good 
agreement between kinematics measured using pPXF and 
GANDALF; the largest offsets are observed in velocity disper- 



sion, and are only significant in galaxies with strong emission 
lines and relatively weak continuum flux. 

GANDALF also allows for spectra to be cleaned of de- 
tected emission. Here we follow Sarzi et al. (2006) and adopt 
an amplitude-to-noise (A/N) threshold of 4, where noise is 
deflned as the scatter about the best-fit absorption template. 
As both [Nl] and H/3 emission strongly affect our Lick index 
measurements the removal of these lines results in a signifi- 
cant improvement in the quality of stellar population fits for 
emission- line galaxies (see discussion in Section |4|. The ap- 
proximate detection threshold of emission lines in our data 
can be calculated as a function of S/N using our adopted 
A/N and typical line broadening. If we assume a typical in- 
trinsic line width of ~30-50 kms"'^, then our sensitivity is 
~0.6 A for AAOmega data and ~0.9 A for 6dFGS data at 
the median S/N of our galaxy data. 

We use the Ha/H/3 Balmer decrement to calculate the 
extinction correction for our emission-line fluxes where ap- 
plicable using the Rv = 3.1 reddening curve of Cardelli, 
Clayton & Mathis (1989) and an intrinsic Hq/H/3 line ratio 
of 2.85 (Osterbrock 1989). In galaxies with very low extinc- 
tion the errors on our line-flux measurements can sometimes 
result in a negative extinction estimate. As this is clearly 
non-physical we adopt E{B — V) < 0.01 as an upper limit of 
extinction and assign this to galaxies with low or spurious 
(i.e. negative) extinction measurements. 



3 LICK INDEX MEASUREMENTS AND 

CALIBRATION TO THE LICK/IDS SYSTEM 

Here we use Lick absorption features (Burstein et al. 1984; 
Trager et al. 1998) to characterise the stellar populations of 
our galaxies. Lick indices are defined based on their partic- 
ular sensitivity to either age or metallicity, and so are useful 
for describing these parameters in integrated spectra. Lick 
line-strength measurements were carried out using the in- 
dex definitions of Trager et al. (1998), supplemented with 
higher-order Balmer line definitions from Worthey & Ot- 
taviani (1997). We omit from our Lick index measurements 
Fe5782, NaD, TiOi and Ti02 as these fall in the wavelength 
range of the AAOmega dichroic where sensitivity is decreas- 
ing rapidly. 

3.1 Velocity dispersion corrections 

In order to properly compare measured indices with SSP 
models it is important that the total broadening of any 
given absorption feature matches the wavelength-dependent 
broadening of the Lick/IDS system. Galaxies with combined 
instrumental and velocity-dispersion broadening less than 
the target Lick resolution are broadened as necessary us- 
ing a wavelength- dependent Gaussian kernel to match the 
description of the Lick/IDS resolution given by Worthey & 
Ottaviani (1997). 

In cases where the combined instrument and velocity 
dispersion broadening is greater than the required Lick/IDS 
system resolution we calculate line-width corrections based 
on iteratively broadened template stars. These line-width 
adjustments are then applied to correct measured indices 
back to the desired Lick/IDS resolution. As a test of this 
method. Proctor & Sansom (2002) compare stellar and 
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Table 1. Lick index corrections for AAOmega data and their 
associated errors. 



Index 


Offset 


^rms 




0.148 


0.426 


USf 


0.053 


0.208 


CNi 


0.018 


0.008 


CN2 


0.024 


0.010 


Ca4227 


0.153 


0.098 


G4300 


-0.029 


0.110 


H7A 


0.244 


0.255 


H7f 


0.099 


0.118 


Fe4383 


0.362 


0.093 


Ca4455 


0.386 


0.049 


Fe4531 


0.025 


0.107 


C4668 


-0.569 


0.197 


H/3 


-0.133 


0.039 


Fe5015 


0.222 


0.074 


Mgi 


0.028 


0.002 


Mg2 


0.036 


0.002 


Mgfe 


0.006 


0.051 


Fe5270 


-0.096 


0.125 


Fe5335 


-0.004 


0.055 


Fe5406 


-0.014 


0.039 


Fe5709 


0.183 


0.090 



galaxy spectra broadened using this technique, finding good 
agreement between the two. For further details see Proctor 
& Sansom (2002). 



3.2 Lick/IDS system calibration 

In order to calibrate index measurements to the Lick/IDS 
system, we have observed several stars from the Lick stellar 
library at the beginning of each night. Index measurements 
for AAOmega data are then adjusted based on a comparison 
of published Lick standard star indices and our stellar ob- 
servations. The mean offsets found for each index are listed 
in Table [T] along with their associated r.m.s. error. 

We are unable to calibrate 6dFGS data in the same way 
as AAOmega data, and so we generate two sets of indices 
for 6dFGS galaxies: the first using the corrections provided 
in Table [TJ and the second with no correction. We refer the 
reader to i]4.1l for a comparison of the properties derived 
using these two sets of indices for overlapping 6dFGS and 
AAOmega galaxies. 



3.3 Index errors 

The fitting technique that we adopt for these data involves 
the simultaneous minimisation for a large number of Lick 
indices (~20). In order for these fits to be reliable, an accu- 
rate representation of index errors is crucial, but notoriously 
difficult to obtain for spectral observations. 

Proctor et al. (2008; hereafter P-l-08) have recently used 
spectra from the 6dFGS DRl (Jones et al. 2004) to measure 
ages and metallicities in a large sample of galaxies using Lick 
indices. In estimating their index errors, P-l-08 have used 
duplicate observations of the same galaxy (from overlapping 
survey regions) to calculate the r.m.s. scatter in their Lick in- 
dex measurements at a given S/N. The r.m.s. value was then 
fit as a function of S/N using the form r.m.s. = a/(S/N-f 6), 




20 40 

S/N per A 

Figure 1. H/3 index error vs. signal-to-noise for galaxies in our 
spectroscopic sample. The solid line shows the best fit to our 
measured index errors as a function of S/N (see i|3.3| l. while the 
dashed line shows a similar fit for 6dFGS DRl data derived by 
P-l-08. The dotted, vertical line shows our adopted S/N cut of 12. 



where a and h are constants, which was then used to calcu- 
late representative errors for the rest of the galaxy sample 
based on their S/N. This method not only encapsulates red- 
shift and velocity dispersion errors, but also errors resulting 
from sky subtraction, fibre flat-fielding and throughput cal- 
ibration that are nearly impossible to quantify in individual 
fibre observations. 

Here we adopt a similar strategy to P-f 08 and exploit 
the composite nature of our final spectra by re-measuring in- 
dices on individual observations, of which there are 12 and 34 
observations each for bright and faint galaxies respectively. 
We then adopt the arms for these repeat measurements as 
an accurate representation of the stochastic errors in our 
sample (i.e. Poisson noise, sky-subtraction etc). In Fig. [T] we 
show how the H/3 index error derived in this manner varies 
as a function of S/N for the galaxies in our sample, and for 
comparison we also show the line describing index errors as 
derived by P-l-08 for 6dFGS DRl data. Error curves for all 
indices are shown in Appendix 1X1 

For our AAOmega data the total adopted index error 
is a combination of this random error, the (Tr.m.s. associated 
with conversion to the Lick/IDS system ( H3.2p and the index 
broadening due to errors on velocity dispersion and recession 
velocity measurements ( ^2.3p . For 6dFGS data we adopt the 
recession velocity and velocity dispersion errors as discussed 
in i)2.3l as well as the index errors as a function of S/N 
derived by P-l-08. 



4 STELLAR POPULATION MODELS AND 
FITTING 

Having compiled relevant Lick index measurements and 
their corresponding error, we now turn to the task of in- 
terpreting these measurements in terms of the stellar pop- 
ulation parameters they describe. Here we use the Single 
Stellar Population (SSP) models of Thomas, Maraston & 
Bender (2003) with account for variable abundance ratios as 
calculated by Korn, Maraston & Thomas (2005; hereafter we 
refer to the Thomas et al. 2003 and Korn et al. 2005 models 
collectively as TMB03). These models cover the metallicity 
range -2.25<[Z/H]<0.65 with ages from 1 to 15 Gyrs. The 
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variable abundance ratios included by TMB03 consider N, 
O, Mg, Ca, Na, Si and Ti as enhanced elements and, while 
these fall outside a strict definition of a-elements, we will 
hereafter refer to enhancement measures from the TMB03 
models as a-element abundances or [a/Fe]. 

It is well known that spectra are degenerate with re- 
spect to age and metallicity in so much as old, metal-rich 
and young, metal-poor stellar populations are remarkably 
similar. In order to minimise the influence of this degener- 
acy on our final measurements, we adopt the minimisa- 
tion technique of Proctor et al. (2004), simultaneously fit- 
ting as many Lick indices as possible. This method has been 
shown to recover reliable ages, metallicities and a-element 
abundances in both galaxies (e.g. Proctor & Sansom 2002; 
Brough et al. 2007; Proctor et al. 2008) and globular clusters 
(Beasley et al. 2004; Pierce et al. 2006; Mendel et al. 2007). 

Our goal is to obtain a stable fit between our galaxy 
data and the TMB03 models, and to accomplish this we 
adopt an iterative approach to our fitting. As a first step 
we perform a fit using as many indices as our data allow, 
generally between 18 and 21. In subsequent iterations indices 
deviant from the best-fit solution at the 5, 4 and 3 cr levels 
are clipped. A final, manual inspection of the fit is then 
carried out to ensure that stable fits have been attained and 
that indices have not been over- or under-clipped as a result 
of poor error estimates. 

Figs. [2] and [3] show the mean deviation of indices from 
the best fit model for AAOmega and 6dFGS data respec- 
tively. Galaxies free of emission are generally fit very well by 
the TMB03 models and we clip at most 2 indices from any 
given galaxy. Emission-line galaxies are significantly more 
difficult to fit, however the emission-cleaned spectra output 
from GANDALF result in a large improvement to the over- 
all quality of fits for AAOmega data (Fig [2l3,c) and slight 
improvement for 6dFGS data (Fig. [3}3,c). 

Final errors on our stellar population measurements are 
estimated using a Monte Carlo method which re-samples 
the best-fit model convolved with our observed index errors. 
These error estimates represent our best attempt to quantify 
the random error contribution to our measurements, how- 
ever there is an additional systematic error present in the 
models themselves which we do not account for (see, e.g., 
Conroy, Gunn & White 2008). 

The 6dFGS and AAOmega spectra used in this work 
can only be used to derive luminosity-weighted ages, metal- 
licities and a-element abundances. This is important to 
any subsequent analyses as these luminosity-weighted val- 
ues are particularly affected by small centralised bursts of 
star formation and may not be representative of the global 
galaxy properties. In addition, because of aperture size dif- 
ferences between AAOmega and 6dF (2.1" vs. 6.7"; 0.3 kpc 
vs. 0.9 kpc at the distance of NGC 5044) we expect some 
variation in the derived stellar population parameters due 
to both sampling differences and the presence of radial gra- 
dients, predominantly in metallicity. Final measurements of 
stellar population parameters (ages, metallicities etc.) and 
their associated errors for the 67 galaxies in our spectral 
sample are included in Appendix IB] 
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Figure 2. Mean deviation of indices from the best fit TMB03 
SSP model. Open squares show the mean deviation for all indices 
in each galaxy, while filled circles show the mean deviation only 
for indices included in the fits, (a) Galaxies with no detectable 
emission, (b) galaxies with visible emission in their spectra and (c) 
emission galaxies with spectra cleaned using the GANDALF routine 
(see ^Zit . 
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Figure 3. Same as Fig. [2] but for 6dFGS data. Data shown are 
without correction to the Lick/IDS system (see i|4.1| l. 



4.1 Comparison of AAOmega and 6dFGS 

galaxies: agreement and aperture effects 

In Fig. m we show the comparison of stellar population 
parameters derived for overlapping AAOmega and 6dFGS 
galaxies. In this comparison we have used indices measured 
from 6dFGS data with no Lick/IDS system correction ap- 
plied. Using indices adjusted using the AAOmega Lick/IDS 
system correction gives similar results, but with an increase 
in scatter. We therefore use uncorrected 6dFGS indices for 
the remainder of this work. 

Of particular interest here is the comparison between 
galaxy properties sampled using the different fibre aper- 
tures of AAOmega and 6dF. Perhaps the most obvious trend 
is observed in metallicity, which is also the most robust of 
our stellar population measurements. 6dFGS galaxies show 
a systematic offset towards lower metallicities of ~0.1dex, 
consistent with observations that the majority of galax- 
ies possess negative gradients in metallicity (i.e. galaxies 
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Figure 4. Comparison of metallicity ([Z/H]), a-element abundance and age for tlie 6 overlapping 6dFGS and AAOmega galaxies with 
S/N > 12. Open (red) and filled (black) circles represent galaxies with and without emission lines. Dashed lines in each panel represent 
equality. Parameters for 6dFGS data are shown with no applied correction to the Lick/IDS system (see i|3.2ll . 



are more metal-poor at larger radii, e.g. Sanchez-Blazquez 
et al. 2007; Brough et al. 2007). 

Comparisons of a-element abundances and ages are 
considerably more scattered than metallicity, indicative of 
the increased difficulty in their measurement. Broadly, we 
expect the larger fibre aperture of 6dF to yield older ages, 
particularly for galaxies with recent or currently ongoing 
star formation as these bursts are (generally) centrally con- 
centrated. We find this to be the case for half of the 6dFGS 
galaxies, but the other half scatter to younger ages with 
considerable uncertainty. As these galaxies are just above 
our adopted S/N cut, we consider these measurements to be 
highly uncertain rather than representative of a real trend 
in the data. We expect the variation in aperture size be- 
tween AAOmega and 6dF to have little systematic effect on 
the measured a-element abundances, as observed gradients 
in [a/Fe] are both weak and variable. In Fig. |4l a-element 
measurements using AAOmega data are consistently higher 
than those from the 6dFGS data. Again, this is most likely 
a S/N effect as the small dynamic range of a-element abun- 
dances makes their determination difficult in low S/N data. 

Motivated by the comparisons in Fig. |4] and the dis- 
cussion above, in forthcoming sections we consider galaxy 
metallicities to be robust for both the 6dFGS and AAOmega 
samples. Age and a-element abundances for 6dFGS galaxies 
can be used to give a rough indicator of stellar population 
trends at larger group-centric radii (see i|6.1l and Paper I), 
but we will refrain from using them to make any specific 
judgements due to the large uncertainty in their relation- 
ship to AAOmega data. 



5 GALAXY PROPERTIES 

5.1 Stellar population parameters versus stellar 
mass 

Before attempting to place our observed galaxy population 
in the context of the NGC 5044 group environment, we first 
examine the properties of galaxies independent of their par- 
ticular place in the group. In Fig. [S] we show how our mea- 
sured ages, metallicities and a-element abundance vary with 
stellar mass (M,). Stellar masses have been computed using 
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Figure 5. Relationship between stellar population properties and 
stellar mass (see il5.H l. Filled (black) and open (red) symbols rep- 
resent passive and emission-line galaxies (see i|2.4l l. with circles 
and triangles delineating galaxies measured using AAOmega and 
6dFGS spectra, respectively. Errors bars represent 1-sigma devi- 
ations estimated using the Monte Carlo technique described in 
Section |4] 
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GALEXEV stellar population synthesis models of Bruzual 
& Chariot (2003; hereafter BC03). Here we use BC03 mod- 
els constructed using the "Padova 1994" isochrones (see 
BC03 and references therein) and stellar initial mass func- 
tion (IMF) of Chabrier (2003). Using our derived central 
ages and metallicities, we then extract the appropriate B- 
band mass-to- light ratio (M/L) from the BC03 models. 

Of the three stellar population parameters measured 
and shown in Fig. [S] metallicity shows the strongest cor- 
relation with mass, spanning more than two decades with 
a scatter of ~0.26 dex. The general tightness of the mass- 
metallicity relation as determined from both emission- and 
absorption- line analyses is well known from large sam- 
ples such as the SDSS (e.g. Tremonti et al. 2004; Gallazzi 
et al. 2006), however our data probe this relation an order 
of magnitude lower in both mass and metallicity. We see 
evidence for increased scatter among emission-line galax- 
ies ('^ 0.39 dex for emission-line galaxies vs. ~ 0.23 dex for 
passive galaxies), which could be indicative of "contami- 
nation" from recently formed stellar populations, however 
the fact that metallicities for low-mass emission-line galaxies 
(log Mt < 9) appear depressed relative to their passive coun- 
terparts suggest that contamination is not playing a signif- 
icant role, as recently formed populations should be metal- 
enriched. Instead, the observed low metallicities (as well as 
increased scatter) could be an effect of metal-enriched out- 
flows in these galaxies, driven either by starbursts or super- 
novae, which could prevent the incorporation of metals into 
newly-formed stars. We note, however, that measurements 
for emission-line systems are also affected by the quality of 
our applied emission-line corrections, which would naturally 
lead to an increased scatter in absorption-line measurements 
in these systems. 

An interesting question related to the mass-metallicity 
relation is the degree to which it depends on environment. 
While a detailed analysis of this dependence requires a 
comparison of galaxies at fixed stellar mass in a range of 
environments, even a qualitative comparison of the mass- 
metallicity relation in several different environments should 
allow identification of any gross systematic deviations. To 
this end, we have selected a comparative sample of galax- 
ies, both at higher and lower masses, in two other environ- 
ments, the Local Group and the Shapley supercluster. Data 
for Local Group dwarf galaxies are taken from the compila- 
tion of Woo, Courteau & Dekel (2008), who provide stellar 
masses, derived using the colour- M,/L relation of Bell & 
de Jong (2001), and average [Z/H] from red giant branch 
stars. Metallicities for a sample of galaxies in the Shapley 
supercluster have been calculated using Lick line index mea- 
surements from Smith et al. (2007) and the fitting methods 
described in Section|4l Stellar masses for these galaxies have 
been estimated as described above, using B-band magni- 
tudes from Smith et al. (2007). A comparison of the mass- 
metallicity relation for these three samples is shown in Fig. 
[6l In this figure, metallicity estimates for Shapley superclus- 
ter galaxies have had additional offset of -1-0.15 dex applied 
to account for the different physical aperture size between 
NGC5044 and Shapley observations (0.3 kpc vs. 1.9 kpc re- 



spectively )0. Data shown in Fig. [6] form a remarkably uni- 
form mass-metallicity sequence over roughly six orders of 
magnitude in galaxy stellar mass, despite the widely varied 
environments sampled. Furthermore, apart from the obvious 
differences in stellar mass range, the data suggest that en- 
vironment is likely playing only a small role in establishing 
the mass-metallicity relation, consistent with recent results 
from the SDSS (e.g. Mouchine et al. 2007; van den Bosch 
et al. 2007). 

Interpretations of the mass-metallicity relation gener- 
ally focus around decreasing star-formation efficiency at 
low masses, which can suitably explain the trends ob- 
served in galaxies (e.g. Tremonti et al. 2004; Savaglio 
et al. 2005; Gallazzi et al. 2006) and is reproduced in hydro- 
dynamic simulations including feedback mechanisms (e.g. 
Brooks et al. 2007). In this interpretation of the mass- 
metallicity relation we expect to see low-mass galaxies trend 
towards lower a-element abundances, as the relationship of 
Q-element abundance to supernovae timescales, namely SNII 
versus SNIa contributions, provides leverage in separating 
between rapid, high-efficiency star formation and on-going, 
low-efficiency star formation (e.g. Terlevich & Forbes 2002; 
Thomas et al. 2005). 

Examining the observed relation of [a/Fe] with mass in 
Fig. O we find considerable scatter, particularly in galax- 
ies with emission. The majority of this scatter results from 
increasingly large index errors at low S/N (see Appendix 
\K\ . but there is also a systematic effect related to our 
fitting of low-mass, and hence low-metallicity, data. The 
TMB03 model grids "pinch" together at low [Z/H], making 
the discrimination of [a/Fe] increasingly difficult (Mendel 
et al. 2007). So, while we appear to observe significant 
scatter in the [a/Fe] measurements of low-mass galaxies, 
which would generally indicate significant variation in star- 
formation histories, we are unable to draw any strong con- 
clusions from these data. If we restrict ourselves to galaxies 
with log M* > 8.5 Mq, then our data exhibit a weak positive 
correlation of a-element enhancement with mass at the ~95 
percent level, consistent with the interpretation of the mass- 
metallicity relation as indicative of star-formation efficiency, 
discussed above. 

In our age estimates we see a clear offset between 
emission-line and passive galaxies, unlike either metallicity 
or a-element abundance. The mean stellar age of emission- 
line galaxies is ~3.6 Gyrs, relative to ~9.1Gyrs for passive 
galaxies. We must, of course, interpret these age measure- 
ments with caution as they are strictly luminosity-weighted 
measurements within our fibre apertures. Emission-line 
galaxies show a trend of younger central ages towards in- 
creasing mass, significant at the ~2(t level, which we believe 
to be largely driven by aperture effects and the centrally 
concentrated nature of star formation. Analyses of central 
stellar populations in conjunction with global photometric 
measurements allow one to constrain the relative fraction of 
mass contained in central starbursts, which has been found 
to be of order 10 percent for a range of galaxy masses (e.g. 
P-(-08). If we consider the fact that the fraction of galaxy 
light sampled by AAOmega's 2" aperture varies by more 



^ Assuming an average metallicity gradient of -0.25 dex per dex, 
e.g. Sanchez-Blazquez et al. (2007) 
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Figure 6. Mass-metallicity relation for galaxies spanning a broad range of environments. Filled circles, open circles and open squares 
represent data from the NGC 5044 group, Local Group and Shapley supercluster respectively. Errors for the Local Group data are taken 
as the average error quoted in Woo et al. (2008) of 0.17 and 0.2 dex for stellar mass and metallicity. Errors on NGC 5044 group and 
Shapley supercluster metallicities have been estimated using the Monte Carlo technique described in Section |4] 



than an order of magnitude from the brightest to faintest 
sources in our sample, then the stellar population measure- 
ments for our most massive galaxies are almost completely 
dominated by any central star formation, while in smaller 
galaxies we sample a growing fraction of the underlying, old 
stellar populations. 

The passive sample of galaxies seems to separate into 
two sub-samples: one with relatively uniform, old ages 
( > 9 Gyrs) , and a second population with relatively young 
central ages. This younger sub-sample almost certainly rep- 
resents galaxies which have more recently undergone bursts 
of star formation and are now fading to older apparent ages. 
One possible cause for this population could be recent, gas- 
rich mergers which would serve to drive down central age 
measurements (e.g. Kauffmann 1996). In this scenario we 
would expect to see these star-formation bursts accompanied 
by a decrease in a-element abundance as new generations of 
stars form from increasingly metal-enriched gas, however the 
a-element abundances for these galaxies are consistent with 
the bulk of the passive galaxy population. 

To summarise, the ages, metallicities and a-element 
abundances of galaxies in the NGC 5044 group are consis- 
tent with properties found in larger galaxy samples. The 
tight mass-metallicity relation and its interpretation as a 
sequence of star-formation efficiency is consistent with our 
measurements of a-element abundance, although the trend 
of elemental abundance with mass appears relatively weak 
in these data. 



5.2 The age— mass— metallicity relation 

The comparison of galaxy ages, metallicities and element 
abundances described in the previous section demonstrates 
the close relationship between metallicity and stellar mass, 
but it should also be considered that more complex relation- 
ships between multiple stellar population parameters and 
mass may exist. For example, in a study of local early- 
type galaxies Trager et al. (2000) found that their ellipti- 



cal galaxies were described by two separate two-dimensional 
relations in four-dimensional space: a plane described by 
the linear combination of logt, log a and [Z/H], and a re- 
lationship between [a/Fe] and log cr such that more mas- 
sive galaxies have higher values of a-element enhancement. 
Smith et al. (2008b) found a similar relation between age, 
mass and metallicity for dwarf galaxies in the Coma cluster, 
noting in particular an increase in scatter at low masses. 

Motivated by recent comparisons suggesting a tight re- 
lationship between the properties of galaxies at a fixed stel- 
lar mass (e.g. van den Bosch et al. 2008), we fit the plane 
described by our data for age, mass and metallicity such 
that 

log M, = a [Z/H] -f /3 log t -f 7, (1) 

using a least-squares fitting method and minimising resid- 
uals orthogonal to the plane. The best-fitting parameters 
for these data result in a plane such that a = 1.83 ± 0.46, 
(3 = 1.41 ± 0.79 and 7 = 9.00 ± 0.61. In Fig. [7]we show pro- 
jections of this best-fit plane in stellar mass, M*, metallicity, 
[Z/H], and age. There is a clear separation between passive 
and emission-line galaxies about the fitted plane, which is 
driven primarily by the difference in mean age of the two 
populations and evident in Fig. [T}:. The separation between 
these two populations is inconsistent with being an effect of 
the age-metallicity degeneracy in our fits, which primarily 
moves galaxies along the trends in Fig. [T] If we consider fits 
to the emission- and non-emission-line galaxies separately, 
both populations show a similar dependence on stellar mass; 
the offset between these two populations is dominated by a 
varying age-metallicity relation. Low-mass galaxies in our 
sample exhibit an increase in scatter about the fitted plane 
relative to their high-mass counterparts, in agreement with 
the observations of Smith et al. (2008b). 

The projections shown in Fig. [7] are helpful to examine 
the distribution of galaxy parameters relative to one an- 
other, but are somewhat difficult to interpret in physical 
terms. In Fig. [8] we show a more standard projection of the 
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Figure 8. Age— metallicity relation for NGC 5044 group galaxies. 
Symbols are the same as in Fig. |5] Lines of constant stellar mass 
are shown for logM» = 8, 9, 10, 11 Mq. Error bars are the 1-sigma 
errors determined from Monte Carlo simulations (see 

plane in terms of age and metallicity, where dashed lines 
represent the age-met allicity relation at fixed stellar mass 
as derived from the best-fit plane discussed above. 

If we naively interpret the distribution of galaxy ages 
and metallicities shown in Fig. [S] then there is is little ev- 
idence for a relationship between central galaxy age and 
metallicity; galaxies formed at a similar time in the early 
universe exhibit a spread in metallicity of two or more dex. 
However, our fits in the age-mass-met allicity plane suggest 
that this is a flawed conclusion. While it is true to say that 
there is no clear age-metallicity trend when considering our 
total sample of galaxies, subsamples of data in narrow mass 
bins (one or two orders of magnitude) suggest that galaxies 
with older central ages are more metal poor than centrally 
young galaxies, supporting evolution of the mass-metallicity 
relation over time (e.g. Kobulnicky et al. 2003; Brooks 
et al. 2007; Lamareille et al. 2007; Maiohno et al. 2008). 

5.3 Predicted vs. observed colours 

While up to this point we have used the BC03 models to de- 
rive M/L, and hence mass, for our galaxies, it is important 
to establish the extent to which these SSP models accu- 
rately describe the properties of our sample galaxies. To un- 
dertake this comparison we have supplemented the B- and 
K-hand photometry used in Paper I with B — V and g — i 
colours for a subsample of NGC 5044 group dwarf galaxies 
presented by Cellone (1999) and Cellone & Buzzoni (2001; 
2005; S. Cellone, private communication). Predicted colours 
are calculated using the BC03 models and measured ages 
and metallicities. Comparisons of these predictions to our 
observed colours are shown in Fig. 

Galaxies lacking in emission show relatively good agree- 
ment with the predicted BC03 colour for their age and 



metallicity. This serves as an excellent confirmation in the 
reliability of M/L estimates for these galaxies which, to some 
extent, can be seen from the tightness of the mass-metallicity 
relation in Fig. [5] The scatter in colours of emission line 
galaxies is considerable, and the causes of this are likely to 
be two-fold: firstly, we are sampling fundamentally different 
regions with our spectroscopy and photometry. Whereas this 
has a negligible effect for passive galaxies, the centrally con- 
centrated star-formation in our emission galaxies means that 
the colours we predict should be, and indeed are, generally 
bluer than the observed "global" photometry. In addition, 
stellar population parameters derived for our emission-line 
galaxies are the most uncertain, which leads to a greater un- 
certainty in both our measurements of ages and metallicities 
and hence a corresponding uncertainty in their predicted 
colours. Overall, however, the above comparisons suggest 
that our galaxy data are reasonably well described by the 
BC03 models. 



6 GLOBAL PROPERTIES AND GALAXY 
DISTRIBUTION 

We have so far examined general trends in the galaxy popu- 
lation of the NGC 5044 group in their own right. We now 
turn to a more general discussion of the distribution of 
galaxy properties within the group. 

6.1 Stellar population trends with radius 

By combining semi-analytic models for star formation and 
galaxy evolution with recent, large-scale Af-body simulations 
such as the Millennium Simulation (Springel et al. 2005), 
recent theoretical work has made great strides in predict- 
ing the distribution of galaxy properties in massive struc- 
tures like groups and clusters. As an example, de Lucia 
et al. (2006) found that the luminosity-weighted age, metal- 
licity and stellar mass fall with increasing distance from 
the cluster centre in their models, out to the virial radius. 
Given the hierarchical formation scenario favoured in cur- 
rent ACDM models these results are not surprising; those 
galaxies in the highest density regions form first, and subse- 
quent generations of galaxies accreted by the cluster are dis- 
tributed with radius according to the redshift at which they 
become cluster members (e.g. Gao et al. 2004). It should 
be noted, however, that there is considerable scope for the 
interpretation of projected cluster-centric distance as an in- 
dicator of accretion time to be confused by three-body in- 
teractions which can eject bona fide members to beyond the 
virial radius (see e.g. Ludlow et al. 2008). 

In light of these predictions, we are motivated to ex- 
amine the radial distribution of stellar populations in the 
NGC 5044 group. Galaxies residing in groups and clusters for 
more than a dynamical time should undergo some degree of 
mass segregation due to the increasing efficiency of dynam- 
ical friction with total galaxy mass. In apparently relaxed 
groups and clusters, such as NGC 5044, we expect any such 
segregation to be readily apparent. In Paper I we examined 
the distribution of dwarf and giant galaxies in the NGC 5044 
group using _B-band magnitudes as a simple discriminator 
between low- and high-mass systems but found no signifi- 
cant evidence for differing radial distributions of these two 
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Figure 7. Projections of the fitted plane to age, mass and metaUicity. Symbols are the same as Fig. \5\ Dashed lines represent equality 
and are shown for reference. 
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Figure 9. BC03 predicted colours vs. observed global B — K, B — V and g — i. Filled and open symbols represent passive and emission-line 
galaxies respectively. Errors on predicted colours are calculated using the errors estimates on our central ages and metallicities. Dashed 
lines represent equality between the predicted and observed colours. 



sub-populations. Here, we are able to refine this analysis 
using the stellar masses derived from the BC03 models. 

In the upper panels of Fig. [TD]we show the distribution 
of stellar mass, age, a-element abundance and metaUicity in 
fixed radial bins. In the lower panels we show the radial dis- 
tribution of galaxies in bins of stellar mass, age, a-element 
abundance and metaUicity. This figure includes only those 
galaxies within 500 kpc of the group centre (roughly two- 
thirds the group's virial radius) as this is approximately the 
region observed uniformly with our new AAOmega obser- 
vations, and thus less prone to spurious, selection-induced 
trends. 

Focusing first on the distribution of properties in fixed 
radial bins, we find no significant differences in the mean 
values for age or a-element abundance (shown by vertical 
dashed lines in the top panels of Fig. llOp . Data show some 
evidence for a deviation in mean galaxy stellar mass and 
metaUicity, and a Kolmogorv-Smirnov (KS) test confirms 
that metallicities in the inner 150 kpc of the group are in- 
consistent with being drawn from the same distribution as 
either of the outer two radial bins at the 2.5 sigma level. 
Despite the apparent offset in mean stellar mass, a KS test 
does not find galaxies in the inner group region to have a sig- 
nificantly different mass distribution relative to other radial 
bins. 



Turning to the radial distribution of galaxy proper- 
ties (Fig. 1101 lower panels), a KS test finds that both 
low-mass and low-metallicity galaxies (logM* < 8.8 and 
[Z/H]< — 0.7dex) exhibit radial distributions that are differ- 
ent to their higher mass and metaUicity counterparts. Qual- 
itatively, these differences are evident in Fig. [TO] by the clear 
lack of galaxies from the lowest mass and metaUicity bins in 
the inner 100 kpc of the group centre. The fact that we see 
this segregation in both mass and metaUicity is expected 
from the tight correlation between these two quantities in 
our data; the converse holds true for the lack of peculiar 
distributions in age and a-element distributions given their 
lack of a strong correlation with mass. 

In order to investigate the relationship between mass, 
metaUicity and radius further, in Fig. [11] we show galaxy 
metallicities plotted against their group-centric radii, where 
galaxies have been separated into three separate mass bins. 
This figure serves as a visual confirmation of the mean offset 
of galaxy metaUicity observed within ~ 150 kpc in Fig. 1101 
and also shows that this trend is due primarily to a lack of 
low-mass, and hence low-metallicity, galaxies from the inner 
100 kpc of the group. Finally, Fig. [11] suggests that there are 
no strong radial trends in metaUicity at fixed stellar mass. 

Previous work by Mathews et al. (2004; hereafter 
M+04) using the NGC 5044 group catalogue of Ferguson 
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Figure 10. The radial dependence of various properties of galaxies in the NGC 5044 group. Top panels show the distribution of stellar 
mass, age, «-element abundance and metallicity in three different radial bins of ij < 150 kpc, 150 kpc < R < 250 kpc and R > 250 kpc 
(solid orange, green and blue hatched regions respectively). Vertical dashed lines show the mean of each bin. Lower panels show the 
radial distribution of galaxy properties in three separate bins of stellar mass, age, o-element abundance and metallicity, where bins have 
been selected to contain roughly equal numbers of galaxies. 



& Sandage (1990; hereafter FS90) has noted a similar lack 
of low-luminosity galaxies in the inner group region. As- 
suming dwarf galaxies act as tracer particles in the group 
potential, M-l-04 compared the cumulative number distribu- 
tion of dwarf galaxies to the predicted surface mass distribu- 
tion of a Navarro-Frenk- White (NFW) dark matter profile 
(Navarro, Frenk & White 1995, 1996, 1997), finding that 
the NGC 5044 group falls well below the NFW prediction 
inside ~300 kpc. M-l-04, and later Faltenbacher & Mathews 
(2005; hereafter FM05), offer two explanations for this ap- 
parent central deficit of dwarf galaxies. M-)-04 argue that a 
lack of dwarf galaxies could be tied to a low survival rate 
of low-mass galaxies entering the NGC 5044 group halo at 
high redshift (2 — 3 < 2 < 6), or due to suppression of star 
formation in low-mass galaxies as a result of AGN or out- 
flow activity from the central galaxy. The follow-up work by 
FM05 showed, using a simple dynamical model, that tidal 
disruption and dynamical friction could also be responsible 
for the apparent lack of dwarfs at small radii. In their model, 
low-mass galaxies approaching the central group potential 
could be strongly disrupted or destroyed. 

Our sample has allowed us to refine the NGC 5044 group 
membership used by M-(-04, and we find that limiting the 
FS90 group samples to only the spectroscopically defined 
members in Paper I does not significantly alter the conclu- 
sions of M-|-04j; there still appears to be a deficit of dwarf 
galaxies at low projected radii (< 300 kpc). M-l-04 suggest 
that, as a result of low dwarf galaxy survival rates at high 
redshifts, there should be a relative lack of old, low-mass 
galaxies in the group. The data presented in Fig. [10] do not 
suggest any such trend in galaxy ages, with low-mass galax- 
ies appearing similar in mean age to their high-mass counter- 
parts (see also Figs.[5]and[8]). In addition, the M-l-04 scenario 
would suggest that, relative to massive galaxies, low-mass 

Where the NFW profile is scaled to match the cumulative 
galaxy distribution at ~350 kpc as in M-l-04 



dwarf galaxies should be kinematically "younger" in their 
velocity distribution in the group; in Paper I we find no 
such evidence for kinematic segregation between high- and 
low-luminosity galaxies. A repeat analysis of sub-population 
kinematics using the stellar masses derived in this work gives 
a similar result, further suggesting that the scenario posed 
by M-l-04 does not describe these data in full. 

Perhaps the most plausible explanation is that put for- 
ward by FM05, where the apparent central deficit of dwarf 
galaxies is a result of both dynamical friction and tidal dis- 
ruption. In their work, FM05 find that the evolution of 
galaxies with stellar to total mass ratios of ~20 best de- 
scribe the population of disrupted galaxies in the NGC 5044 
group. Using the relationship between stellar mass and total 
baryonic mass determined from fits to star-forming galaxies 
in the SDSS (Baldry, Glazebrook & Driver 2008), we can use 
the cosmic baryon fraction, = 0.171 as determined from 
the WMAP 5yr results (Komatsu et al. 2008) to estimate 
the total halo mass for galaxies of a given stellar mass. For 
the stellar-to-total mass ratio predicted by FM05 of 20, the 
above above relation gives a corresponding stellar mass of 
log A/* ~ 8.6. This is in agreement with the lower limit of 
stellar mass that we observe in the central 100 kpc of the 
group centre, providing support for the tidal disruption sce- 
nario of FM05. 

Confirmation of this would need to come from addi- 
tional deep imaging to observe the low surface-brightness 
intragroup light which, given the possibly large number of 
disrupted or destroyed dwarf galaxies, could be as much as 
~35 percent of the group's total luminosity (FM05). 

6.2 Galaxies and the hot intragroup medium 

The metal content of the intracluster/intragroup medium 
(ICM/IGM) is directly related to the star formation his- 
tories of the galaxies via outflows and supernovae driven 
winds, and therefore provides a useful diagnostic of both 
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Figure 11. Galaxy mctallicity against projected group-centric 
radius. Galaxies are divided into three mass bins with approx- 
imately equal numbers such that circles, squares and triangles 
represent galaxies with logAf* > 9.6, 8.8 < logM* < 9.6 and 
logM, < 8.8. Solid, dashed and dot-dashed lines represent the 
mean metallicity for each mass bin. 



group and galaxy evolution. Despite the presence of ongo- 
ing evolution in the stellar populations of group and cluster 
galaxies, observations of high redshift X-ray samples suggest 
that the bulk of ICM metals were already in place at 2; ~ 1 
(e.g. Tozzi et al. 2003; Balestra et al. 2007) and therefore ne- 
cessitate that most star formation activity and subsequent 
ICM enrichment takes place at early times. 

Rasmussen & Ponman (2007; hereafter RP07) have used 
Chandra archival data to analyse the X-ray abundance prop- 
erties of 15 nearby groups, including NGC 5044. From their 
data RP07 extract radial profiles of iron and silicon abun- 
dance in the hot IGM out to ~ 200 kpc in NGC 5044, facil- 
itating a comparison between our galaxy stellar population 
measurements and the chemical properties of the intragroup 
gas. In Fig. [12] we plot the X-ray derived radial iron and sili- 
con abundances along with the stellar mass-density profiles. 

In their abundance analysis of hot gas in the NGC 5044 
group, RP07 show that the central silicon to iron ratio is 
roughly solar (i.e. [Si/Fe]~ 0.0), suggesting SNIa are play- 
ing an important role in enriching gas in the central group 
regions. While we lack significant overlap between the X- 
ray profiles and our binned group profile, the slopes in both 
iron and silicon abundance in the intragroup gas are consis- 
tent with that of the stellar mass density between 50 and 
~100kpc. In clusters, the central "excess" of iron in the 
central regions relative the rest of the ICM is normally at- 
tributed to enriched outflows from the central galaxy (e.g. 
Bohringer et al. 2004) . In the case of NGC 5044 this remains 
a plausible explanation for the high central iron abundance, 
particularly given recent evidence for large outflows asso- 
ciated with NGC5044's central AGN (Temi, Brighenti & 
Mathews 2007; Gastaldello et al. 2007). We note, however, 
that our data are also consistent with satellite galaxies hav- 
ing contributed signiflcantly towards building up the central 
iron peak in the group, particularly given the extended na- 
ture of the central excess, ~50kpc, relative to the optical 
extent of NGC 5044 (~20kpc; FS90; Paturel et al. 2000; 
shown with black dashed line in Fig. I12p and the the good 
correspondence between the radial profile of iron abundance 
and stellar mass-density out to beyond ~ 100 kpc. 

In the outer regions, R > 100 kpc, there is evidence 
for a rising silicon abundance, in contrast to the iron abun- 



dance which continues to fall. While the significance of this 
upturn is only marginal due to its dependence on the outer- 
most measurement in the X-ray data, we have no reason to 
believe this is a spurious measurement. The observed rise in 
[Si/Fe] is in agreement with X-ray observations for numer- 
ous groups and clusters (e.g. Finoguenov et al. 2000; RP07) 
and suggests an increased contribution of SNII relative to 
SNIa at large projected radii. 

To explore the relation between galaxies and the intra- 
group medium further, in Fig. [T3]we show the same binned 
X-ray abundance profiles as in Fig. 1121 only this time over- 
layed on the radial galaxy trends for metallicity and [a/Fe] 
shown in Fig. 1111 This figure clearly shows that both the 
IGM iron abundance and the [Si/Fe] ratio are well below 
the average values measured in galajcies at similar radii, de- 
spite both iron and silicon declining with radius at a similar 
rate to the stellar mass-density. On the other hand, low- 
mass galaxies generally exhibit stellar metallicities which 
fall below the measured IGM iron abundance at a given 
radius. Given the shallow potential wells of low-mass galax- 
ies, metal-enhanced starburst or supernovae driven winds 
provide a plausible explanation for the observed offset (see 
also Section [5. l|l . Low-mass galaxies not withstanding, the 
observed low chemical abundance of the IGM is consistent 
with the findings of Buote, Brighenti & Mathews (2004), 
who also note the low abundance of iron in the IGM rela- 
tive to the group luminosity (compared to clusters), as well 
as other work suggesting a high fraction of "primordial" gas 
in the IGM (e.g. Gibson & Matteucci 1997; Moretti, Porti- 
nari & Chiosi 2003). 

In terms of silicon abundance, we find a surprising cor- 
relation between the rise in [Si/Fe] at ~100 kpc and the pres- 
ence of low-mass, dwarf galaxies at similar radii. If there is a 
connection between these two observations, it would imply 
a very high contribution of SNII relative to SNIa in dwarf 
galaxies, as well as a significant contribution from these low- 
mass systems to the total gas mass of the group. Previous 
studies, however, have shown that, while dwarf galaxies cer- 
tainly contribute to the reservoir of intragroup gas, it is most 
likely not significant, constituting at most 15% of the total 
gas mass in clusters (Gibson & Matteucci 1997). This sug- 
gests that the connection between these observations is likely 
only coincidental, particularly given the relatively low IGM 
iron abundance outside the group core, 10% solar, and the 
likely sub-solar outfiows from dwarf galaxies. 



7 EMISSION LINE PROPERTIES 

Using the sub-sample of our NGC 5044 group with measur- 
able emission lines, we now discuss the emission-line charac- 
teristics of our group galaxies, including star-formation rates 
and nebular metallicities. 

First, in order to separate galaxies undergoing star 
formation activity from those with other strong ionising 
sources, i.e. AGN, we use a standard emission line diagnostic 
comparing the [OIII]a5007/H/3 and [NII]A6584/HQf flux ratios 
(Baldwin, Phillips & Terlevich 1981), shown in Fig. [14] We 
use two predictions to characterise the emission line flux 
ratios of our galaxies. The first is the theoretical maximum 
starburst limit calculated by Kewley et al. (2001; dashed line 
in Fig. I14|) . while the second is the semi-empirical limit of 
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Figure 12. Comparison of stellar mass density {solid black line) 
and X-ray derived iron and silicon abundances (red dashed and 
dot-dashed lines respectively) as a function of radius. Data have 
had a fixed offset applied to be shown on the same figure. X- 
ray data have been binned for clarity. The vertical dashed line 
represents the radial extent of NGC 5044's light profile of ~20 kpc. 
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Figure 13. Comparison of X-ray derived silicon and iron abun- 
dance profiles {dashed lines in each panel) with the distribution 
of galaxy a-element abundances and metallicities. Galaxy sym- 
bols are the same as in Fig. If II X-ray data have been binned for 
clarity. 



pure star formation defined by Kauffmann et al. (2003) us- 
ing a large sample of SDSS emission line galaxies {solid line 
in Fig. I14|) . These lines separate our sample into galaxies 
which are purely star-forming, strongly AGN or ionised by 
a composite source (some contribution from both AGN and 
star-formation), shown with circles, triangles and squares in 
Fig.lMl 

NGC 5044 is known to host AGN activity (Rickes, Pas- 
toriza & Bonatto 2004, Brough et al. 2007) and is identified 
as such in Fig. 1141 but we also identify two other group 
galaxies as hosts of potential AGN or other strong ionising 
sources: NGC 5037 (FS 068) and FS 082. The classification 
of these galaxies as AGN sources is supported by the avail- 
able Chandra data, in which all three galaxies are detected. 
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Figure 14. Line diagnostic for emission galaxies in our NGC 5044 
group sample. The solid and dashed lines represent the semi- 
empirical pure star formation line from Kauffmann et al. (2003) 
and the theoretical maximum starburst line from Kewley 
et al. (2001) respectively. Symbols represent galaxies which are 
likely purely star-forming, AGN or composite AGN and star- 
forming (circles, open triangles and open squares respectively). 



7.1 Star formation rate 

In large samples, comparisons of galaxy properties with lo- 
cal environment show strong evidence for a decreasing frac- 
tion of star-forming galaxies in regions of higher projected 
galaxy surface density (Lewis et al. 2002; Gomez et al. 2003; 
Poggianti et al. 2006). Even at relatively low gas densities, 
such as those found in groups, recent work has shown that 
ram-pressure stripping can influence star-formation via de- 
pletion of a galaxy's hot-gas reservoir (SivakofT et al. 2004; 
Machacek et al. 2005; Rasmussen et al. 2006; Kawata & 
Mulchaey 2008). However in higher-density environments, 
such as the Coma cluster, Poggianti et al. (2004) have found 
evidence for young post-starburst galaxies preferentially lo- 
cated near the edges of X-ray substructures. These observa- 
tions suggest that galaxy-ICM interactions may play a key 
role in the truncation of star formation, but may equally 
contribute to initiating bursts of star formation in galaxies 
as they encounter dense environments. 

In our own data, varying fibre throughput and sensi- 
tivity make proper fiux calibration of our data difficult and 
therefore limit our ability to derive absolute measurements 
of the star formation rate, SFR, using Ha line fiuxes. In- 
stead, we use the relation derived by Guzman et al. (1997) 
to estimate star formation rates using B-band luminosity 
and [OH] equivalent widths, where 



SFR(M0yr" 



2.47 X 10"^^ Lb(L0)EW 



[Oil] ■ 



(2) 



We use the B-band luminosities for galaxies as published in 
Paper I, which are taken from FS90 and Paturel et al. (2005) . 

The use of [OH] instead of Ha introduces an additional 
error into the estimation of the star formation rate of or- 
der 0.2 to 0.3 dex (Kewley et al. 2004). The relation fit 
by Guzman et al. (1997) to convert between [OH] fiux and 



[on] equivalent width carries with it increased uncertainty 
of ~0.1dex in the derived SFR. Given these two relatively 
significant error contributions, in addition to the relatively 
large errors in the FS90 and Paturel et al. photometry (0.5 
and 0.3 dex respectively, see Paper I), we consider our star 
formation rates to only be proportional to the star forma- 
tion activity in any given galaxy, but by no means a clear 
indicator of the absolute star formation rate. 

In Fig. [TS] we show the specific star formation rate for 
galaxies as estimated using Eqn. [2] in relation to their pro- 
jected radius from the group centre. The star formation rates 
shown in Fig. 1 151 have been normalised by the galaxy stellar 
mass and are shown in units of yr~^ per 10^ Mq. We see no 
clear correlation of SFR with radius, and results are similar if 
we consider projected galaxy density. Motivated by the pos- 
sible dependence of SFR on galaxy-ICM/IGM interaction, 
we have also examined the star formation rate relative to 
a simple expression proportional to the IGM ram pressure, 
Pram = ^(1") | I'group — , whcrc p{r) IS the IGM density and 
It'group — i^l is the galaxy velocity u relative to systemic group 
velocity i^group. Assuming that p(r) oc r"'^, we then arrive at 
a comparison of SFR against jzJgroup — v'^ /r^ , which shows 
a similar lack of correlation as the radial and galaxy den- 
sity comparisons described above. We note that projection 
effects are likely having a strong effect on any of our radial 
SFR comparisons as we cannot properly relate our galaxy 
positions and recession velocities to their true position and 
relative motion withing the group. Nevertheless, while these 
data do not allow us rule out the role of interactions with 
intragroup gas as inffuencing the observed star formation 
in these galaxies, they seem to hint that these effects are 
relatively limited in this case. 

An alternative possibility for triggering or enhancing 
star formation is through galaxy-galaxy interactions. Hy- 
drodynamic simulations modelling gas have shown that tidal 
interactions can enhance or induce central star formation by 
driving gas from the disk to the central regions (e.g. Barnes 
& Hernquist 1992; Di Matteo et al. 2007), which supports 
observations of enhanced star formation in close pairs of 
galaxies (e.g. Lin et al. 2007; Ellison et al. 2008). An exami- 
nation of DSS images for the 22 galaxies with [Oil] emission 
shows that only 2 galaxies, FS90 134 and FS90 137, ex- 
hibit clear evidence for ongoing interaction. Several other 
galaxies show some evidence of disturbed morphology (i.e. 
boxy bulges and disturbed disks), however it appears that 
the majority of recent star formation activity we see in the 
NGC 5044 group is not driven by currently ongoing galaxy- 
galaxy interactions. 

7.2 Emission-line metallicity 

Nebular metallicity and underlying, stellar metallicity are 
intimately related via star formation. The evolution of stellar 
metallicity (as measured by absorption features) is strictly 
dependent on a galaxy's star formation history, while neb- 
ular metallicities continue to vary as previous generations 
of stars evolve and enrich the interstellar medium. We are 
therefore able to use the relative offset of emission- and 
absorption- line metallicities to extract information related 
to the last significant star formation episode. Here we use 
a reparameterisation of the Kewley & Dopita (2002; here- 
after KD02) _R23 method presented by Kobulnicky & Kew- 
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Figure 15. Specific star formation rate as it varies with projected 
distance from the group centre. Errors are indicative of the uncer- 
tainty in our EW[on] ^rid S-band luminosity measurements, but 
do not include the systematic contribution from using [Oil] over 
Ha and equivalent widths as opposed to line fluxes (see i|7.1l l. 

ley (2004; hereafter KK04) to derive gas-phase metallicities 
for a subsample of our emission-line galaxies. We exclude 
from this analysis galaxies exceeding the maximum starburst 
limit shown in Fig. [14] as likely hosting AGN (NGC 5044, 
NGC 5037 and FS082). 

The i?23 metallicity calibrator is typically defined us- 
ing emission line fiuxes, however this requires spectra to be 
properly flux calibrated as line ratios spanning a significant 
wavelength range will otherwise be incorrect. An alternate 
approach is to use measurements of equivalent width in place 
of line fiuxes to compute EW--R23, which avoids the issues 
of flux calibration and extinction correction, but introduces 
the problem of potential variation in the underlying stellar 
continua (particularly between old and young spectra). In 
a statistical sense, the variation between i?23 and EW--R23 
results in an an additional error on the final metallicity de- 
termination of ~0.11 dex (Kobulnicky & Phillips 2003), how- 
ever in our relatively small sample where we know the under- 
lying stellar populations vary significantly, the uncertainty 
between _R23 and EW--R23 can have a potentially significant 
effect on our results. Considering both the wide range of 
spectra types and uncertain ffux calibration of our data, we 
estimate nebular metallicities using both _R23 and EW--R23. 
Since the effects of varying stellar continua and ffux cali- 
bration affect emission-line metallicity measurements differ- 
ently, by using both of these measurements we can assess 
the reliability of our conclusions. 

Comparison of EW--R23 and -R23 shows that the scat- 
ter between the two measures is ~0.08dex, or ~0.11dex in 
terms of the final determination of 12 -I- log(0/H) (Kob- 
ulnicky & Phillips 2003), which is generally less than the 
intrinsic error in the calibration of ii23 to nebular oxygen 
abundance. In our own data, the scatter between EW--R23 
and -R23 is ~ 0.16 dex, but this scatter is primarily driven 
by two galaxies with significant offsets between EW--R23 and 
J?23. Excluding these galaxies the scatter on our data falls to 
~0.07dex. We therefore exclude these deviant galaxies from 
the remainder of our emission line analyses and consider 
measurements on the remaining galaxies to be reliable. 

Comparison of stellar and nebular metallicities is car- 
ried out by converting nebular oxygen abundances to the 
implied metal mass fraction using the relation Z ~ 29 x 
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]^q[i2+(o/h)]-12^ ^ggjjjjjjjjg ^jjg standard solar abundance dis- 
tribution and solar oxygen abundance of 12 + (O/H) = 8.72 
(AUende Prieto, Lambert & Asplund 2001; Kobulnicky & 
Kewley 2004). In Fig. [16] we show the comparison between 
the converted nebular metallicities and stellar metallicity. 
We find that, in general, derived nebular metallicities are 
super-solar, and in some cases nearly 1 dex higher than their 
stellar counterparts. The relatively shallow relation between 
stellar and nebular metallicities is somewhat puzzling, as 
it implies a significantly weaker mass-metallicity relation 
among star-forming galaxies than is generally observed (e.g. 
Tremonti et al. 2004). 

Given this somewhat curious result we are motivated to 
examine the relative difference between nebular and stellar 
metallicity as a function of other galaxy properties. In Fig. 
1171 the relative off'set between nebular and stellar metallic- 
ity is plotted against galaxy stellar mass, age and a-element 
abundance. The relationship between metallicity offset and 
stellar mass in Fig. [17] is driven primarily by the mass- 
metallicity relation and shows that the most massive galax- 
ies have the lowest measured nebular-stellar metallicity off- 
sets as a result of the relatively uniform, high nebular metal- 
licities measured. 

In terms of age, those galaxies with the lowest metallic- 
ity offsets also have the youngest luminosity- weighted ages, 
which is consistent with the observed stellar population 
forming recently from gas present in these galaxies. Inter- 
estingly, while the trend of metallicity offset with age is 
consistent with our understanding of galaxy formation, the 
magnitude of this ofi'set is large, nearly 1.5 dex in the case 
of the lowest-mass galaxies. This offset is difficult to explain 
through pure passive evolution of the stellar populations in 
these galaxies, i.e. through enrichment from SNII and SNIa, 
which in extreme cases can account for ~ldex of metallic- 
ity increase (see, e.g. Sansom & Proctor 1998). An alternate 
explanation is that ongoing star formation has helped to en- 
hance the observed gaseous metallicities through the evolu- 
tion of metal-enriched stars. In this case the integrated stel- 
lar populations should also appear both younger and more 
metal rich in those galaxies hosting the largest nebular- 
stellar metallicity offsets, in conflict with the trends in Fig. 
1171 Finally, these data could be explained in a scenario where 
the gas reservoirs of low-mass galaxies are being preferen- 
tially depleted over high-mass galaxies, which would increase 
the effect of any subsequent enrichment due to supernovae. 

In order to explain these observations in terms of gas 
removal, we require that the removal process must be more 
efficient in low-mass galaxies, ruling out AGN activity as 
a possibility due to its strong scaling with mass. Starburst 
outflows, primarily driven by SNII, satisfy the above stated 
mass dependence but would also likely affect the observed 
[a/Fe] ratio of stars formed from any remaining gas, in 
conflict with Fig. 117b . Galaxy-IGM interactions such as 
ram-pressure stripping (Gunn & Gott 1982) and viscous 
stripping (Nulsen 1982) satisfy both the inverse mass de- 
pendence and non-preferential gas removal (i.e. gas-phase 
metals from SNII and SNIa must be removed more-or-less 
equally) stipulated by these observations. Strictly speak- 
ing, we cannot rule out any processes from contributing to 
the observed galaxy processes, however our data are best 
explained by galaxy-IGM dominating gas removal in the 
NGC 5044 group. 



in 




Figure 16. Comparison of metallicity derived from the nebular 
oxygen abundance, [Z/H]r23 ^nd stellar metallicity, [Z/H]Lick, 
derived using Lick indices. The dashed line represents equality 
between the nebular and stellar metallicities. 

8 DISCUSSION AND CONCLUSIONS 

In this second of two papers, we have undertaken a spectro- 
scopic investigation of galaxies in the group environment, de- 
riving stellar masses, ages, metallicities and a-element abun- 
dance ratios for 67 of the 111 spectroscopically confirmed 
NGC 5044 group members described by Mendel et al. (2008) . 
These measurements have allowed us to examine the star 
formation histories of these group members, both as a pop- 
ulation in their own right and in the context of the group 
environment and their location in it. 

The mass-metallicity relation plays a fundamental role 
in describing galaxy populations, and we see strong evidence 
of this in the overall tightness of the mass-metallicity rela- 
tion in our group data over multiple orders of magnitude 
in stellar mass. In the context of environment, by compar- 
ing the mass-metallicity relation of NGC 5044 group galax- 
ies with galaxies in the local group and the Shapley super- 
cluster, we show that galaxies appear to form a continuous 
mass-metallicity relation across several orders of magnitude 
in system mass and spanning upwards of six orders of mag- 
nitude in galaxy stellar mass. These data argue for a relative 
independence of the mass-metallicity relation from environ- 
ment (e.g. Maiolino 2008; van den Bosch et al. 2008). 

In terms of galaxy ages and a-element abundances, we 
find no evidence for strong trends with galaxy mass. Our 
a-element abundance data do not rule out the interpreta- 
tion of star formation efficiency as the primary driver of 
the mass-metallicity relation. Particularly at higher stellar 
masses (logM, >9M0) and metallicities, where wide sep- 
aration of model grids allow for finer abundance measure- 
ments, our data show hints of a trend of increasing a-element 
ratio with mass, consistent with literature results (e.g. Ter- 
levich & Forbes 2002; Thomas et al. 2005). We have shown 
that the surface described simultaneously by age, mass and 
metallicity is vital to analyses of galaxy data spanning a 
broad range in mass. Most notably, mass dependence is vi- 
tal to interpretation of the galaxy age-metallicity relation. 

In terms of the galaxy radial distribution, we find no 
apparent radial dependence in either age or a-element abun- 
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Figure 17. OflFset between nebular and stellar metallicity plotted against stellar mass, age and a-element abundance. 



dances. A KS test finds that low-mass and low-metallicity 
galaxies liave radial distributions inconsistent with their 
high-mass and high-metallicity counterparts, and an exami- 
nation of this apparent difference shows that we observe no 
galaxies with log M* < 8.9 Mq within '^lOO kpc projected ra- 
dius from the group centre. In the context of previous anal- 
yses of the NGC 5044 group, M-l-04 and FM05 have both 
noted the apparent lack of dwarf galaxies in the group's cen- 
tral regions using the photometric catalogue of FS90. Our 
spectroscopic observations confirm that the dwarf galaxy 
distribution observed by M+04 and FM05 is not due to 
the inclusion of significant foreground or background con- 
tamination in the FS90 catalogue (see Paper I). Of the pos- 
sible scenarios for explaining this deficit, the most plausi- 
ble seems to be that dwarf galaxies are being tidally dis- 
rupted or (icstrojrcd through interactions with the group or 
central galaxy potential. Our data axe consistent with the 
simulations of FM05 in their prediction of stellar-to-total 
mass ratios for disrupted galaxies of ~ 20, corresponding 
to ^ log M, ~ 8.6. Deep widc-fiold imaging is required 
to confirm this conclusion, as the models of FM05 predict 
stars from the disrupted population of dwarf galaxies to con- 
tribute nearly 35 percent of the total group luminosity, con- 
sistent with estimates of intraclustcr light for other groups 
at a similar mass (e.g. Gonzalez, Zaritsky & Zabludoff 2007). 

With regard to similar types of analyses on other groups 
and clusters. Coma is the only cluster analysed to a sim- 
ilar depth with stellar population measurements. Smith 
et al. (2008a) have examined radial gradients of Coma clus- 
ter dwarf galaxies, finding strong trends towards younger 
ages and higher metallicities with increasing cluster-centric 
radii. While these results are at odds with our own, we be- 
lieve this disagreement to be largely driven by the galaxy 
selection employed by Smith et al. (2008a). The trends ob- 
served by Smith et al. (2008a) are consistent with the in- 
terpretation that the south-west "infall" region of Coma 
has already undergone a passage through the cluster centre. 
Observations of younger ages and higher metallicities in a 
"quenched" dwarf population would then be consistent with 
the idea of induced star formation from galaxy encounters 
with the dense/hot IGM (e.g. Poggianti et al. 2004). The 
fact that we do not see such a trend would seem to suggest 
that we either fail to measure galaxies at sufficient radii to 



sec effects of induced star formation in our galaxies, or that 
such dramatic interactions are not occurring within the cen- 
tral regions we probe here. It is worth noting that in Paper 
I we identified a possible infalling substructure to the north- 
east of NGC 5044 at ~1.4Mpc projected separation, but our 
stellar population measurements do not extend that far and 
so we are unable to comment on possible trends at these 
large radii. 

While we do not see any obvious trends of stellar popu- 
lations with radius, as Smith et al. do, we do perhaps find a 
more subtle interpretation of galaxy evolution in the group 
environment through analysis of emission-line properties. 
Comparison of emission- and absorption-line metallicities 
has allowed us to extract information regarding the previ- 
ous star-formation episodes of galaxies through the relation 
of nebular and stellar metallicities. What we find is that 
the relative offset between these two metallicity estimates 
is too large to be explained by either pure passive evolu- 
tion of galaxies, which can account for at most ~ 1 dex of 
offset, or ongoing star formation, which is inconsistent with 
other observed galaxy properties. The mass-dependence of 
our emission-absorption metallicity offset instead is consis- 
tent with a scenario in which galaxies are being stripped of 
their gas primarily through galaxy-IGM interactions, and 
therefore subsequent ISM enrichment is having a larger ef- 
fect on the observed emission-line metallicity. We require a 
larger sample of galaxies to confirm such an analysis, how- 
ever this method provides a novel technique for probing en- 
vironmental effects on galaxy evolution. 

We turn, finally, to the question of NGC 5044's evolu- 
tion as a group and what information the galaxy population 
holds in this regard. The relative lack of radial trends in the 
stellar population data limit the extent to which the gradual 
accretion scheme present in semi-analytic models of group 
formation (e.g. dc Lucia et al. 2006; Berricr et al. 2008) can 
be acting in the NGC 5044 group, at least without consider- 
able mixing of the galaxy distribution with radius (e.g. Lud- 
low et al. 2008). The large peculiar velocity of NGC 5044 it- 
self relative to the group mean (^ 200 kms~^), possible pres- 
ence of substructure on the outer regions of the group (see 
Paper I) and recent detection of a "sloshing" type cold front 
in the group's X-ray distribution (Gastaldello et al. 2008) 
suggest that the group is still growing and in the process 
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of relaxation. However, the group's Gaussian line-of-sight 
velocity distribution and relatively uniform X-ray profiles 
argue against any recent, dynamically violent evolution of 
the group. 

The work presented here represents the first targeted 
effort at constraining the stellar population trends of group 
galaxies. Further work is clearly needed, particularly on 
larger samples of homogeneously selected galaxies, to better 
constrain the processes driving galaxy evolution in groups. 
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APPENDIX A: INDEX ERRORS 

Here we show the error curves for all Lick/IDS indices as 
derived from AAOmega observations. 



20 Mendel et al. 




S/N per A S/N per A 

Figure Al. Index Error as a function of signal to noise. Solid Figure A2. Same as Fig. lAll 

(black) and dashed (red) lines represent the best fit hyperbolic 
functions for these data and the 6dFGS DRl data of P+08 re- 
spectively. The vertical dotted line represents our adopted S/N 
cut of 12. 
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APPENDIX B: STELLAR POPULATION 
MEASUREMENTS 

Here we include the (Iciived stellar mass, M/L, age, metal- 
licity and a-element abundance for all 67 galaxies in our 
spectral sample. 
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Table Bl. Stellar population parameters of NGC 5044 group galaxies. Only those galaxies with S/N > 12 are included (see 

text.) 
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FS90 107 




13 


15 59.30 


-16 


23 


49.80 


10.3+0-'l 


3.73lo;^» 


0.93±0.05 


0.20±0.05 


0.27±0.05 


FS90 108 




13 


15 59.90 


-16 


10 


32.00 


8.6+0-1 

— U. i 


'-7Q~\~0 

— U.oo 


1.07±0.09 


-0.85±0.05 


0.12+0.14 


FS90 117 




13 


16 23.11 


-16 


8 11.35 


9.4+0-0 


9 Oq + 0.17 

•^•""J—O 15 


0.95±0.05 


-0.30d=0.05 


0.18±0.06 


FS90 123 




13 


16 35.62 


-16 


26 5.70 


8.5+0+ 


3.45 Q 22 


1.15±0.06 


-0.68±0.05 


0.09±0.10 


FS90 127 




13 


16 39.40 


-16 


57 3.00 


8.4+0-'l 

— U.I 


3-41+S-ll 


1.02±0.07 


-0.70±0.06 


-0.30±0.10 


FS90 133 




13 


16 55.43 


-16 


26 32.20 


7 7+0.1 
'•'—01 


o.82lo,:?o 


0.53±0.09 


-0.93±0.08 


0.46±0.17 


FS90 134 




13 


16 56.23 


-16 


35 34.70 


9.il?,:l 


n 7r. + 0.08 


0.47±0.07 


-1.25±0.08 


0.06±0.13 


FS90 135 




13 


16 54.40 


-16 


3 18.00 


8.210:1 


2.4710:- 


1.07±0.09 


-0.85±0.06 


0.42+0.23 


FS90 137 




13 


16 58.49 


-16 


38 


5.50 


10.5lO;l 


0.6810:1" 


0.12±0.07 


0.10±0.06 


0.00±0.06 


FS90 138 




13 


17 3.10 


-16 


35 


39.00 


8.010-1 


2 04+0-1'^ 


1.02±0.10 


-1.02±0.05 


0.46±0.23 


FS90 144 




13 


17 10.79 


-16 


13 


46.74 




49+0.23 
o.'±»_j, 22 


1.05±0.05 


-0.17±0.05 


0.30±0.05 


FS90 153 




13 


17 36.38 


-16 


32 


25.40 


9 5+0-0 
^•''-O.O 


2 60+0-14 


0.88±0.05 


-0.33±0.05 


0.06±0.08 


FS90 158 




13 


18 9.04 


-16 


58 


15.60 


8 9+0-1 


2 92+0-19 


1.10±0.06 


-0.97±0.05 


O.OOiO.ll 


FS90 161 




13 


18 17.80 


-16 


38 


25.00 


8 0+0-2 
°-"-o.i 


1 91+0-77 

-^•^^-0.40 


1.02±0.20 


-1.20±0.07 


0.42±0.37 


2MASXJ13102493-1655578 


13 


10 24.96 


-16 


55 


57.50 


9 8+0-2 


1 41+0.32 
^•^^-0.29 


0.70±0.15 


-0.57±0.13 


0.32+0.15 


2MASX.J13114576-1915421 


13 


11 45.77 


-19 


15 


42.30 




4 1 1+0.29 
^-^^-0.30 


1.00±0.06 


-0.05±0.05 


0.09±0.07 


2MASXJ13115849-1644541 


13 


11 58.49 


-16 


44 


54.10 


9-6+^:1 


2 98+0-12 


l.OOitO.06 


-0.53±0.05 


0.00±0.12 


2MASXJ13125449-1645571 (FS90 005) 


13 


12 54.50 


-16 


45 


57.00 


10.610:1 




1.02±0.06 


-0.10±0.05 


0.12+0.07 
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Table Bl— Continued 



Galaxy Identifier 



RA 
(J2000) 



DEC 
(J2000) 



log M« 



M/Le 



log age 



[Fe/H] 
(dex) 



[a/Fe] 
(dex) 



(1) 






(2) 




(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


2MASXJ13130549- 


■1628411 


(FS90 009) 


13 


13 


5.49 


-16 


28 41.72 


"•''-0.1 


q Q'3 + 0-28 
O.^0_g 3Q 


0.85±0.08 


-0.07+0.06 


-0.03+0.10 


2MASXJ13143041- 


■1732009 




13 


14 


30.42 


-17 


32 0.90 


g 4+0.1 


1 Oq + 0.18 

^•^■^-0.12 


0.90+0.07 


-1.40+0.05 


-0.30+0.14 


2MASXJ13143485- 


■1629289 


(FS90 050) 


13 


14 


34.86 


-16 


29 28.90 


8 

S.O_Q 2 


41+0.47 
"•^-^--0.13 


0.10+0.29 


-0.57+0.20 


-0.30+0.12 


2MASXJ13150409- 


■1623391 


(FS90 075) 


13 


15 


4.08 


-16 


23 40.20 


9 4+0.1 

"•^-0.1 


q Oq+0.29 

0.i50_Q 27 


0.97+0.07 


-0.35±0.05 


-0.21±0.09 


2MASXJ13153203- 


■1628509 


(FS90 094) 


13 


15 


32.04 


-16 


28 51.10 


Q 1+0.2 
='•^-0.2 


^•D^-0.34 


0.45±0.14 


0.00±0.14 


-0.18±0.12 


2MASXJ13164875- 


1620397 




13 


16 


48.75 


-16 


20 39.70 


9 5+°-i 


2.93+0-2? 


1.07+0.07 


-0.75+0.05 


0.09+0.16 


2MASXJ13165533- 


1756417 




13 


16 


55.35 


-17 


56 42.00 


9 5+0. 1 


^-^-^--0.28 


0.97+0.08 


-0.88+0.05 


-0.06+0.14 


2MASXJ13165624- 


1635347 


(FS90 134) 


13 


16 


56.23 


-16 


35 34.70 


9 6+0-1 


2 50+0-22 


1.07+0.09 


-1.50+0.06 


-0.30+0.18 


2MASXJ13171239- 


■1715162 




13 


17 


12.40 


-17 


15 16.10 


9 7+°l 
"•'-0.1 


1 34+0-18 

^••^^-0.13 


0.62+0.10 


-0.78±0.08 


-0.15+0.11 


2MASXJ13184125- 


■1904476 




13 


18 


41.26 


-19 


4 47.70 


0.9_Q 2 


39+0-18 


-0.05±0.18 


-0.10±0.07 


-0.06±0.15 


2MASXJ13185909- 


■1835167 




13 


18 


59.09 


-18 


35 16.70 


10 i+°-i 


q 91+0.24 
•^•^-^-0.25 


0.85±0.06 


0.03±0.05 


0.09±0.07 



